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Coherence of single spins coupled to a nuclear spin bath of varying density 
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The dynamics of single electron and nuclear spins in a diamond lattice with different ^^C nuclear 
spin concentration is investigated. It is shown that coherent control of up to three individual nuclei 
in a dense nuclear spin cluster is feasible. The free induction decays of nuclear spin Bell states and 
single nuclear coherences among ^^C nuclear spins are compared and analyzed. Reduction of a free 
induction decay time T2 and a coherence time T2 upon increase of nuclear spin concentration has 
been found. For diamond material with depleted concentration of nuclear spin, T2* as long as 30 /is 
and T2 of up to 1.8 ms for the electron spin has been observed. The ^^C concentration dependence 
of T2 is explained by Fermi contact and dipolar interactions with nuclei in the lattice. It has been 
found that T2 decreases approximately as 1/n, where n is ^^C concentration, as expected for an 
electron spin interacting with a nuclear spin bath. 
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Defect centers in diamond have attracted considerable 
interest recently owing to their application for quantum 
information processing, communication and metrology. 
[l|,i,i,i,i,li,l3 Especially the nitrogen- vacancy (NV) 
center, with its strong and spin dependant optical transi- 
tions allows for single spin readout and exquisite coherent 
control which is crucial for quantum information applica- 
tions. [J, [1, 0, 05 [5"] Owing to the high Debye temperature 
of diamond and weak coupling to acoustic phonons NV 
electron spins show long coherence time. It was e.g. pro- 
posed to build small quantum registers by exploiting the 
interaction between the electron spin and a small num- 
ber of nuclear spins in the immediate vicinity. Five-qubit 
would be sufficient to perform all functions necessary for 
a node in a defect center based quantum repeater node. 
[1, m Up to now coherent control, swapping of coherence 
and even entanglement between up to two nuclei and the 
electron spin was demonstrated. [31] To increase the size 
of the quantum register, more nuclei need to be coupled 
to the electron spin. The approach taken here is to in- 
crease the concentration of paramagnetic ^^C nuclei in 
the lattice. We systematically demonstrate coherent con- 
trol of up to three nuclear spins being coupled to an NV 
center electron spin in ^^C isotopically enriched crystals, 
notwithstanding the fact that the electron decoherence 
time T2 linearly scales with the ^^C concentration. Fur- 
thermore, our experiments provide experimental insight 
into long studied problem of single central spin coupled 
to a paramagnetic environment. 8, 9, Owing to possi- 
bility to address individual electron spins in matrix with 
adjustable nuclear spin content we show the transition 
from diluted to dense spin bath (the situation relevant 
for spins in GaAs quantum dots). 



The quantum system used in the present work is the 
negatively charged NV center in diamond, which com- 
prises a substitutional nitrogen atom with an adjacent 
vacancy. (Fig. 1(h)) The electron ground state of it is a 
spin triplet. Upon optical excitation the NV center shows 
strong fluorescence allowing it to be observed on an in- 
dividual basis by confocal microscopy. The fluorescence 
intensity of the defect is spin-dependent owing to spin 
selective relaxation via singlet state, which allows opti- 
cal read out of the single electron spin resonance (ESR) 
PJJ] and an efficient electron and nuclear initialization 
at room temperature. Microwaves (MW) and radiofre- 
quency (RF) fields are used for coherent manipulation 
of single electron and nuclear spins using conventional 
ESR/NMR techniques. All measurement is carried out at 
~20°C. Diamond enables for the unique opportunity to 
control the concentration of paramagnetic nuclear spins. 
The most abundant ^^C has zero nuclear spin. The con- 
centration of ^^C nuclear spins (7=1/2) can be adjusted 
to the suitable value by controlling the isotopic content 
of ^^C in the growth medium. Two types of synthetic 
diamonds were used in this study. Crystals with 0.35%, 
1.1%, 8.4%, 20.7% ^^C concentrations were synthesized 
by a MW plasma-assisted homoepitaxial chemical va- 
por deposition technique (CVD) Ij, [l3| with ^^CH4 or 
^2CH4 gases. The 0.03% ^^C and -100% ^^C enriched 
diamonds were synthesized by high pressure high temper- 
ature method (HPHT) using getters preventing incorpo- 
ration of nitrogen into lattice. In all samples the concen- 
tration of paramagnetic impurity measured by ESR was 
under detection limit (below 1 ppb) except 0.03% and 
100% ^^C diamond where the concentration of nitrogen 
was at 1 ppm level. 
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FIG. 1: (Color online) ESR spectra of single NVs with (a) 
zero, (b) one, (c) two, (d) three ^^C in the 1st shell. Blue 
solid lines are simulation spectra. ESR spectra with (e) one 
and (f) two ^^C in the 3rd shell, and (g) one ^^C in the 1st 
shell and one ^^C in the 3rd shell. Blue solid lines are fitting 
curves with Gaussian shape, (h) Atomic structure of NV 
center. The numbers 1, 2, and 3 mark C in the 1st shell. 

NV centers with different numbers of ^^C atoms in the 
immediate vicinity of the electron spin have been inves- 
tigated in 8.4% ""^^C diamond. Fig. l(a-d) shows single 
ESR spectra indicating the hyperfine coupling (HFC) of 
the electron spin to zero, one, two, and three nuclear 
spins in the first shell (fig. 1(h)). A magnetic field of 
^83 Gauss was oriented along the NV axis ([111] -axis). 
To characterize the spin quantum states associated with 
the transitions in Fig. l(b-d), the spectra were simulated 
by exact diagonalization of the spin Hamiltonian 

H = QePeSB + SDS + {SAh - 9nPnhB) (1) 

Here the electron spin S = 1 and ^^C nuclear spins in 
1st shell are taken into account, is the Bohr and 
j3n the nuclear magneton, respectively. Reported val- 
ues for zero field splitting (ZFS) parameter (\D\ = 2.87 
GHz), [15] isotropic electron and nuclear Zeeman g- values 
(^/e=2.0028, ^n=1.40483), and HFC parameters of A\\ = 
205 MHz and A±= 123 MHz [16] with angle of 106° be- 
tween principal axes of ZFS and HFC, yield precise fits 
of the experimental spectra. The small splittings in the 
central signals of the spectra in Fig. l(c,d) are explained 
by a 2nd order perturbation approach. |3| The smaller 
amplitudes in higher frequency are due to absorption of 
MW by wire on the sample. 

In the ESR spectra of 8.4% ^^C diamond, basically two 
types of couplings are immediately visible (Fig. l(b-g)). 



those around 130 MHz originating from first shell ^^C and 
those around 14 MHz. In Csv symmetry, the number of 
equivalent atoms in close shells around NV is 3 or 6. In 
recent theoretical study, [l^ the 14 MHz splitting are 
assigned to ^^C at 3 and 6 equivalent sites in the 3rd 
shell (see Video in [18].) From measuring more than 250 
individual centers and comparing the probability to find 
the 14 MHz splitting with the one predicted from theory 
we assign this splitting to nuclei in the 3. shell, [l^ 

Individual nuclei in the spin cluster around the electron 
are addressed via their particular NMR frequency. Given 
the increase in spectral density apparent from Fig. 1 one 
might wonder in how far individual nuclei remain ad- 
dressable. However, coherent control even in dense spin 
clusters remains feasible as demonstrated in Fig. 2. Even 
in cases where there are three ^^C in the 1st shell i.e. in 
total four qubits, Rabi nutations of single nuclear spins 
can be driven by an additional RE as shown in Fig. 2(c). 
That is because even multiple RE transition frequencies 
originating from nuclei at equivalent positions split due 
to higher order HFC contributions as shown above. The 
selectivity is not limited to the relatively large splitting in 
the 1st shell but can be applied to ^^C nuclear spins in the 
3rd shell. A spectrum of single NV which has one ^^C in 
the 1st shell and one ^^C in 3rd shell is shown in Fig. 1(g). 
We labeled the four nuclear spin states as |00), |01), |10), 
and 1 11) as shown in Fig. 2(a). Rabi oscillations between 

1 00) and |01) could be observed as shown in Fig. 2(d). To 
obtain a similar Rabi frequency for ^^C nuclear spins in 
the 3rd shell as for those in the 1st shell, about 1 x 10^ 
times higher RF power was necessary. This can be in- 
terpreted mainly by hyperfine enhancement^ [19] which 
predicts that Rabi frequency is almost proportional to 
HFC and the square root of the RF power. A 10 times 
smaller HFC in the 3rd shell supports this interpretation. 
A figure of merit which characterizes the quality of co- 
herent control was derived by swapping quantum states 
among individual nuclei. It was estimated by transfer- 
ring polarization back to a detectable electron spin state 
that 82±5 % of polarization was transferred from |00) to 

1 01) for ^^C nuclear spins in the 3rd shell. 

We are now in a position to engineer simple quantum 
states in the spin cluster around the electrons spin. Bell 
states = 1/V^(|00) ± |11)) and = 1/V2(|01) ± 
1 10)) can be generated from the four effective nuclear spin 
states in Fig. 2(a). In the present case they were prepared 
from two ^^C spins at the 1st shell. The procedure fol- 
lows previous studies [3, |20|] and is schematically shown 
for $~ in Fig. 2(k). After its generation, ^~ shows a 
free induction decay (FID) which is made visible with an 
entanglement detector sequence (Fig. 2(k)). In 8.4 % -"^^C 
diamond, the free induction decay times T| of ^~(T2*^^) 
and <I>~(T2(^^)) were estimated to be 22.0 ± 3.0/is and 
13.3±l.l/is, (see Fig. 2(f,g)) respectively, by fitting with 
exp[-(t/T2*)^] cos(Acc;t), Q where Aco is the detuned 
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FIG. 2: (Color online) (a) Energy levels for nuclear spins in 
M5=-l. Rabi oscillation of single ^^C in the 1st shell of neigh- 
bors around the NV with (b) two ^^C in the 1st shell, (c) three 
^^C in the 1st shell, (d) Rabi oscillation of single ^^C in the 
3rd shell with one ^^C in the 1st shell and one ^^C in the 3rd 
sheh. The pulse sequence is 7r(MW)-Rabi(RF)-7r(MW). [1] 
The ESR transitions of the MW tt pulse are those at lowest 
frequency in Fig. 1 (c,d,g) and are indicated by blue arrows 
in (a). The NMR transitions of the RF pulse are indicated 
by orange arrows in (a). The recording the data in (b-d) re- 
quired about 20 minutes of averaging, (e) ENDOR spectrum 
of ^^C at 1st shell in (b) with the pulse sequence of 7r(MW)- 
7r(RF)-7r(MW). FID of (f) (g) (h) nuclear coherence 
between |00)^|10) and (j) |00)^|01). (k) Pulse sequence for 
$~ generation and detection between two nuclear spins. E 
and Ni,2 marks the electron and the two nuclear spins, re- 
spectively. Spin selective pulses are represented by squares, 
operating on a target qubit. Vertical lines represent logical 
connections. The control qubit state |1) and the state |0) are 
displayed as filled (•) and open (o) circles. For example, (o) 
indicates that the pulse is applied to the target qubit if the 
quantum state of the controlling qubit is |0). 



frequency of FID. As expected from the view point of 
decoherence free subspaces, [21,] a longer T*^^^ compared 
to that of T*^^^ is observed. 

The difference among T2*^^^and T^(^^^ is best analyzed 
when compared with T2* of a nuclear quantum coher- 
ence among states |00)^|10) and |00)^|01). Those 
coherences are labeled as single quantum coherences 
SQl and SQ2, respectively. Their T2* are measured to 



be T. 



2(5Q1) 



=41.1 ± 3.1/is and T*^sQ2) = ^^'^ ^ lAiis, re- 
spectively (Fig. 2(h,j)). The difference of ^2*(5Q1) 
-^2*(5Q2) ™ght be caused by a spatially inhomogeneous 
magnetic noise around the defect caused by an inho- 
mogeneous distribution of ^^C around the two ^^C in 
the 1st shell. Each spin-spin interaction between nu- 
clear spin k surrounding the two ^^C in the 1st shell 
with quantum numbers m/i and m/2 cause oscillation 
given by exp[— i(Aa;i m/i + /\uo2mi2)t\. [22] Here 
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FIG. 3: (Color online) ESR spectra of single NV in (a) 20.7 %, 
(b) 8.4 % ^^C diamond with fitted Gaussian lines (blue), (c) 
Fourier transformed spectrum of FID of 1.1% ^^C diamond 
shown in (d) on expanded frequency axis. The splitting is due 
to HFC of distant ^^C nuclear spins. The hyperfine splitting 
to N is not visible in this spectrum due to polarization of 
the N nuclear spin. [Ssl] (e) Dependence of Inhomogeneous 
linewidth on ^^C concentration. The error bars indicates the 
distributions measured. 



Acji and Acc;2 are spin-spin interaction frequencies of 
the two ^^C in the 1st shell due to surrounding ^^C 
nuclear spins. This implies that ^2*^^^^ and T*^^^ can 

be approximated by l/T^^^) = \l/T^^sQi) - V^2*(5Q2)I 
and l/T^^^)=l/T^^sQi)^^l^2{SQ2Y respectively. Insert- 
ing the measured values for 1/^2*(5Qi 2)' results are 



^2 calc.(^)" 



=25.8ll^/is and T2%^i^^^)=11.4lJ:^/is, respec- 
tively, in good correspondence with measured values. 

Besides the effect of the nuclear spin bath on individ- 
ual ^^C spins the static interaction between the single 
NV electron spin and its environment for different ^^C 
concentrations was investigated by measuring T2* i.e. the 
inhomogeneous ESR linewidth. It was observed that T2* 
increases i.e. the linewidth narrows, with decreasing ^^C 
concentration as shown in Fig. 3. In 0.03 % ^^C diamond, 
an extremely long T2* of 30 /is was found (see Fig. 3). In 
the low ^^C concentration region (^1.1%), the linewidth 
W (full width at half maximum) is derived from T2* by 
W = 2V/n2/7rT2*. The corresponding 18 kHz linewidth is 
the narrowest ever observed for an electron spin in a solid 
material. In the high ^^C concentration region (>1.1%), 
the linewidth is derived from fitting the ESR line of a 
single NV with a Gaussian lineshape. Average values are 
plotted as squares in Fig. 3(e). 

A likely cause for the inhomogeneous ESR linewidth is 
HFC to ^^C nuclear spins. In e.g. Si, the dependence of 
the inhomogeneous linewidth of P donors in ^^Si is well 
fitted by the isotropic HFC {ai) due to Fermi contact 
interaction with ^^Si nuclear spins with a concentration 
(/), W = 2Vm2[fY.i{ail2fY'^. in,!! The sum runs 
over all nuclear spin sites I. In Fig. 3(e), the solid line for 
high ^^C concentrations is calculated by summing only 
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over all the 9 sites in the 3rd shell with ai = 14 MHz 
(see above for assignment of sites and HFC constants). 
It should be noted that contributions from ^^C in the 1st 
shell were not considered in the linewidth calculations 
because they contribute to an observable splitting but 
not to the linewidth. As seen from Fig. 3(e), it fits the 
experimental results well for high ^^C concentration. 

For lower ^^C concentration, experimental data devi- 
ate from this behavior. This is due to the fact that the 
probability that any ^^C is located close to the NV cen- 
ter is getting small upon reduction of ^^C concentration. 
Furthermore, the unpaired electron spin density rapidly 
decreases with distance from the three dangling bonds 
around the vacancy. This is known from the HFC param- 
eters 16|,[l3 which indicates that almost 100 % spin den- 
sity is localized on the C sites in the 1st and the 3rd shell. 
That is why in this situation the most prominent con- 
tribution to the inhomogeneous linewidth is the weaker 
dipole-dipole interaction between electron spin and ^^C 
nuclear spin at distant sites. The lower line in Fig. 3(e) 
is the linewidth 



W = V(/i0Me/in^e^n/4^/^)'(3.195 X lO^^n), (2) 

calculated from the 2nd moment [25] with more than 
3,000 lattice sites for each ^^C concentration (n). Con- 
tributions from ^^C in the 1st and 2nd shell are not con- 
sidered. As seen from Fig. 3(e), W fits the experimental 
results in the low ^^C concentration (^1%) quite well. 
Obviously at low ^^C concentration the linewidth is dom- 
inated by dipole-dipole interaction. 

The dephasing time T2 of the electron spin is measured 
by two pulse Hahn echo decay curves (Fig. 4(a-c)). We 
analyzed T2 of the diamond made by CVD and excluded 
the 0.03 % and 100 % ^^C diamond made by HPHT be- 
cause paramagnetic impurities could not be suppressed 
in HPHT. In the 1.1 % ^^C diamond, a T2 of 0.65 ms was 
found, which is the longest T2 in diamond with natural 
abundance of ^^C measured so far and for the lower ^^C 
concentration of 0.3% an even longer T2 of 1.8 ms was 
measured. [26] T2 is found to be inversely proportional 
to the ^^C concentration as plotted in Fig. 4(d). 

In a theoretical analysis of T2 by the disjoint cluster 
approach, [5] the relationship of T2 ~ (CA(^~^I^ is pro- 
posed, where C is the averaged nuclear-nuclear dipolar 
interaction in the bath and Ac is some characteristic value 
for the electron-nuclear dipolar interaction. Since both 
interactions scale linearly in ^^C concentration (n), T2 
decreases approximately as 1/n in this model. The fitted 
line to the data shown in fig. 4 (d) supports this inverse 
proportionality and fits our data. Note that our data also 
fit the values calculated in [5|] within 30% deviation. 

In conclusion, coherent control of up to three individ- 
ual nuclei in a dense nuclear spin cluster is demonstrated. 
The ^^C concentration dependence of T2* and T2 of elec- 
tron spin point towards ^^C nuclei as the main cause for 
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FIG. 4: (Color online) Echo decays of electron spin in (a) 
20.7 %, (b) 8.4 %, and (c) 1.1% ^^C diamond. The MW pulse 
sequence is 7r/2-r-7r-T-7r/2 where r is delay [8]. Red lines are 
curves fitted with exp[-(t/T2)^]. (d) Plot of T2 over C con- 
centration n. The solid line is fitted with a 1/n dependence. 



dephasing in otherwise clean diamond. The correspon- 
dence with the theoretical line of T2 @ is very important 
to elucidate the dephasing mechanism and to make T2 
longer for quantum information devices [1] and magne- 
tometry. [y, 0] Furthermore, the results show that the 
threshold (^10^ operation) for quantum error-correction 
schemes [27| can be exceeded even in ^^C enriched dia- 
mond at room temperature with typical single-qubit flip 
of several ns. 
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